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1 Introduction 

According to cosmological observations, today our universe has the accelerated expansion 
[1, 2]. What produces that acceleration is the so-called dark energy. Moreover, there 
is dark matter inferred from the existence of gravitational effects on visible matter and 
background radiation, but is undetectable by emitted electromagnetic radiation [3, 4]. Recent 
observations imply that today 72% of all the energy densities of the universe comes from dark 
energy, while dark matter contributes to 23%. The remaining 5% is well-known matter and 
radiation [5, 6]. 

There are many models of dark energy. The most useful and simple model is cosmo- 
logical constant. The results of dark energy from cosmological constant can properly fit 
the observations, but it encounters cosmological constant problem and coincidence problem. 
Since the energy densities of cosmological constant and dark matter are significantly different 
throughout the history of the universe, the problem is posed why their energy densities are of 
the same order at the present time. That problem is referred to the coincidence problem. The 
dark energy models with varying energy density are proposed in order to solve the problems 
of cosmological constant. There are many models explaining dark energy and most of them 
are based on the use of scalar field [7, 8]. However, an uncoupled scalar field cannot account 
for the coincidence problem [9]. One of the possibilities to handle the coincidence problem 
is that the scalar field should couple to dark matter [10, 11] and the coupling can lead to an 
accelerated scaling attractor solution [12]. 

However, the experimental evidences of the fundamental scalar particle have not been 
discovered. Moreover, at a more fundamental level there is no reason to exclude the possibility 
of some other higher form field to be the dark energy. These higher form fields can exist in 
some theories such as the string theory. The presence of them does not necessarily violate 
the cosmological principle. For these reasons an effort has been made in using a vector 
field, a one-form field, to play a role of dark energy [13-16]. However, most of the vector 
field models encounter instabilities [17]. Generalizations to higher form fields have also been 
proposed [18, 19]. Two-form field models also have the same problem as vector field models. 
On the other hand, three- form field models have no such problem: they are stable [19]. 
Therefore, it is of interest to consider a three-form field as a candidate for dark energy. 
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Recently, dark energy from three-form field has been studied and it is found that three-form 
field can provide the accelerating universe [20]. However, the coincidence problem cannot be 
solved since there are no stable fixed points at which energy densities of dark matter and 
dark energy are comparable. 

According to interacting scalar models, there are three classes of coupling form, Qj oc 
<t>Pc, Qn oc H{a ( j ) p,f > + a c p c ), and Qui oc {ot^p^ + a c p c ) where subscript <f> and c denote 
the contribution for the scalar field and cold dark matter respectively. For the first type, it 
cannot provide the standard evolution of the universe while with the same parameters the 
coincidence problem is solved [21]. For the second and third types with exponential potential 
and = case, the coincidence problem cannot be solved since there are no stable fixed 
points at which Q^/Qc ~ 0(1) [12]. However, in the case of ^ and uj^ = constant, 
there are the fixed points which provide the possibility to solve the coincidence problem. 
Nevertheless, they also yield the negative energy density of the dark energy [22, 23]. From 
many investigations of the interacting scalar field dark energy, it suggests that it is not easy 
to solve the coincidence problem with the proper behavior of the evolution. Therefore, it is 
worthwhile to investigate whether the interacting three-form field can provide the possibility 
to solve the coincidence with the proper behavior of the evolution and this is the aim of our 
paper. 

It is important to note that the coupled quintessence with growing matter, a matter 
consisting of particles with an increasing mass, provides one of the possible solutions to 
solve the coincidence problem [24]. The candidate of the growing matter is suggested as 
neutrinos and then the observed dark energy density and the equation of state parameter 
are determined by the neutrinos mass. Significant results of this model are investigated and 
compared with the observational data [25]. 

We review dark energy from non- interacting three- form field in section 2. Our analysis 
suggests that the potential of the three-form field should be runaway. This suggestion is 
obtained by ensuring the well-behaved evolution of the universe. Analogous to scalar field 
dark energy, three types of simple coupling form of the three-form field with dark matter 
are examined for various runaway potentials in section 3. In order to provide the consistent 
perturbations in the models, three types of the coupling form are introduced with the ability 
to express in the covariant form. It is found that the coincidence problem cannot be solved. 
To deal with the coincidence problem, we introduce another covariant coupling form in sec- 
tion 4. It is found that this coupling form can provide the ability to solve the coincidence 
problem with the well-behaved evolution. The results of our investigation are summarized 
and discussed in section 5. 

2 Dark Energy from Three-form Field 

Due to the homogeneous and isotropic universe in the large scale, we consider the flat FRW 
universe with the line element 

ds 2 = -dt 2 + a?(t)dx?. (2.1) 

We assume that the universe consists of matter (both baryon and cold dark matter), radiation, 
and dark energy. We describe each of these species as a fluid. In this paper, the fluid is 
assumed to be perfect, the fluid with no viscosity and momentum density. For the perfect 
fluid, the energy momentum tensor takes the form 

T» v = Diag(-p,p,p,p), (2.2) 
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where p is the energy density of the fluid and p the pressure. The Einstein field equation 
gives two equations 



H 



2 



a) 3 



H = --^{p i+Pi \ (2.4) 

i 

where H is the Hubble parameter and k 2 = 8ttG. The conservation of the energy momentum 
tensor for each species can be expressed as 

Pt + 3H(jn+Pi)=0. (2.5) 
Using (2.3) and (2.4), we obtain the acceleration equation 



\ = -y + 3ft) = -y 1>(1 + 3w<), (2.6) 



where is the equation of state parameter for each species. Therefore, the accelerating 
universe due to each species requires that 

0J % < -1/3, (2.7) 

where, ^matter = and cj ra di a tion = 1/3. Thus both matter and radiation cannot yield the 
accelerating universe that we observe nowadays. In order to explain the current accelerating 
universe, a new constituent of the universe is introduced namely Dark Energy. There are 
many models of dark energy which use scalar field such as quintessence, tachyon, phantom 
[26] and chameleon [27, 28]. Recently, three-form field is introduced as a candidate of dark 
energy models [20]. We will review the dynamical properties of three-form field dark energy 
model in this section. These properties are also compared with the results of the quintessence 
model. We comment on the drawbacks and advantages of the model from the point of view 
of observational compatibility. 

We are interested in a theory of minimal coupling to gravity. The action of a three-form 
field A pup is then 

Sa = - J d*xV^(^R - ^F 2 ~ V(A 2 )) , (2.8) 

where R is the Ricci scalar and V(A 2 ) is the potential of the field and A 2 = A^ up Ap Up . Fp Upa 
is the field strength tensor of the three-form field defined as 

F^ pCT = 4V [/1 ^ H , (2.9) 

where the square bracket means antisymmetrization. Varying the action (2.8) with respect 
to A^p, we get the equation of motion 

V a F a ^ up = YlV\A 2 )A pvp , (2.10) 

where V'(A 2 ) = dV/d (^4 2 )- The Einstein field equation is obtained by varying the action 
(2.8) with respect to g^, 

G^iv = Rfiv — -^Rg^u = K 2 T pu , (2-11) 
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where G^ v is the Einstein tensor, is the Ricci tensor, and is the energy momentum 
tensor given by 



= \F mPl Ff^ + 6V'(A 2 )A^Af - 9/JV 



±F a ^ p F a ^ + V(A 2 )). (2.12) 



For the FRW metric (2.1), where a three-form field depends only on time, a zero component 
of three-form field is nondynamical. From its equation of motion (2.10), it has algebraic 
constraint 12V' (A 2 )Aq pv = 0. Therefore we focus on the space-like components only. In this 
spacetime with the homogeneity and isotropy assumed, the space-like components take the 
form 

A ijk = a 3 (t)e ijk X(t), (2.13) 

where X is a comoving field. The field X is related to the field A fJ,up via A 2 = 6X 2 . From 
(2.10) the equation of motion of the field X is then 



X 



-3HX - V, x -3HX 



(2.14) 



where overdot denotes derivative with respect to the time t and V,x = dV/dX. The energy 
density and pressure of the three-form field are given by 



Px -- 
VX = T\ 



r ° = ^{x + mxf + v{x) 



1 



(X + 3HX) 2 - V(X) + V, x X. 



(2.15) 
(2.16) 



The equation of state parameter of the three-form field is 

Px . . V, X X 

uj x = — = -H • 

Px px 



(2.17) 



Conveniently, we will consider the dimensionless variables for characterizing the dynamics of 
the three-form field. From equation (2.4), (2.5) and (2.14), the dynamical equations are 



3(\/-y-x), 



y' 


= -\Kx)z 2 {xy - W 


!>+ 




4 2 2^ 

-u +v )y 


V 


3 

= — -^v[l + \{x)xz 2 


-w 2 - 


4 2 
~3 U - 


v% 


w' 


3 

= —-w[l + \(x)xz 2 


-w 2 


- A -u 2 - 
3 


v\ 


u 


3 r 4 \ / \ 2 

= + \{x)xz 


-w 2 


4 2 
~3 U ~ 





(2.18) 

(2.19) 
(2.20) 
(2.21) 
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where the dimensionless variables are defined as 



X = KX, y = —( X ' + 3X), 
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-1 


all A 


(c) 




< y* < 1 


y/1 ~ Vl 











-1 


A = 



Table 1. Fixed points in the three- form held dark energy model. 



where pb is energy density of baryon, p r is energy density of radiation, and p c is energy 
density of cold dark matter. The prime denotes the derivative with respect to N = In a, the 
e-folding number. Equation (2.3) is not the dynamical equation. It is the constraint equation 
expressed as 

y 2 + z 2 + u 2 + v 2 + w 2 = 1. (2.24) 

We can find the fixed points of the system by setting x', y' , v', w', and u' equal to zero. We 
summarize all interesting fixed points in Table 1. 

Now we have all equations for characterizing the dynamics of the three-form field. The 
first interesting point is the equation of state parameter of three- form field in equation (2.17). 
It satisfies the condition of accelerating universe in (2.7). Moreover, it can mimic the cos- 
mological constant behavior in the case of vanishing of the potential slope. The equation of 
state parameter can also be adjusted as the phantom phase universe (cj < —1) in order to 
satisfy the recent observational data, u = —1.1 ± 0.14 [29]. Actually, it can take any real 
value. 

Next, we will consider the dynamical equations analytically. From the fact that there 
are matter- and radiation-dominated periods, this yields the result that the three-form field 
dark energy needs to sub-dominate in those periods. This argument will suggest us what the 
shape of the potential is. By considering the matter-dominated period (setting w = 0, y = 
0, u = 0, z = 0), one can get the fixed point giving that period from the dynamical equations 
(2.19) - (2.22). However, from equation (2.18), the value of x is not constrained. At y = 0, it 
grows exponentially when the time is back to the past (N is more negative). This means that 
x must be a huge value at that time. At the matter dominated-period, z needs to vanish. 
This leads to the vanishing of the potential at a large value of X. This suggests us that the 
potential of three-form field should be runaway. However, if we look at the simple runaway 
potential such as an exponential potential, it will not vanish at the finite X. We will discuss 
this issue by taking w = 0,y = 0, u = and z = e where e is a small number. The dynamical 
equations (2.21), and (2.22) vanish automatically. The dynamical equations (2.19) and (2.20) 
can be written as 

y> = y|A(x) e 2 , (2.25) 
v ' = -^l-e 2 (l + xA(x))e 2 . (2.26) 
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In order to ensure that there is the matter-dominated period, y' must be positive and v' must 
be negative at the time after matter-dominated period (If the time is back to the past, it is 
the time before matter-dominated period). For x>l, one obtains the constraint 

A > 0. (2.27) 

This constraint is valid even when we consider radiation-dominated period. Again, the 
constraint of A suggests us that the three-form potential should be runaway. 

We note that the runaway potential we mentioned means that the value of the potential 
at the large X approaches zero lim V{X) — > and the slope of the potential is always 

X— >oo 

negative V^x < 0. The runaway potential will guarantee the well-behaved evolution of the 
background. However, some kinds of non-runaway potential will also provide the well-behaved 
evolution. For example, the power law potential, V = Vox n , with n = 1 also provides the well- 
behaved evolution while the Landau-Ginzburg potential, V = Vq(x 2 — x^) 2 , cannot provides 
such well-behaved evolution. We will consider this issue in detail. Generally, we can set the 
value of z 2 and y 2 to be a small number e 2 and e 2 respectively. Therefore equations (2.25) 
and (2.26) can be generalized as 




v> = -^l-e 2 -e 2 (l + | + xX{x))e% (2.29) 

For the power law potential, V = Vqx u 1 it ensures that the well-behaved evolution will 
be obtained if n < 1 where A = —n/x. Actually n can take the value larger than 1 depending 
on the value of £ 2 /e 2 z - We can see that the non-runaway potential, < n < 1, can also 
provide well-behaved evolution. 

For Landau-Ginzburg potential, V = Vq(x 2 — x^) 2 , with A = — Ax/(x 2 — Xq), one cannot 
get the well-behaved evolution at the large X. Nevertheless, one can argue that the value 
of X will be a small value, (x <C xq), if \j2j?>y always grater than x inferred from equation 
(2.18). However, it is difficult to achieve this situation since y decreases faster than x when 
the time is back to the past during three-form dominated period. We also note that we check 
this analysis with the numerical calculation and found that it agrees with this qualitative 
analysis. 

At the matter- and radiation-dominated periods, we see that the three-form field X is 
a large value. This leads to a huge negative value of X as y oc X + 3HX = 0. Instead a 
scalar field rolls down from the small value like quintessence field, the three-form field rapidly 
climbs up the potential from the large value of field to the small one and reaches the fixed 
point at which dark energy dominates at the present time. One of the most interesting things 
in the three-form dark energy model is the effective potential. From the equation of motion 
of three-form field (2.14) and equation (2.4), the effective potential can be defined through 
the relation 

V eff ,x= (1- ^x 2 )V, x -^x^ liPi . (2.30) 

i 

At the dark energy dominated period, we can neglect the last term which is the contributions 
from matter and radiation. This means that the fixed point at the extremum of the potential 
is x = ^/2/3 which is the same as we consider from dynamical equation (2.18). This is a 
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mechanism by which the potential parameter controls the dynamics of three-form field. It is 
similar to the analysis in reference [20]. Note that in [20], for the potential V = (X 2 - C 2 ) 2 , 
the dynamics of three-form field is controlled by C which is less or more than a/2/3. 

Now, let us analyze the stability of the fixed points. The fixed point (a) corresponds to 
the matter-dominated solution and gives w tot = 0. The eigenvalues are (—3, 3/2, 3); therefore, 
it is not stable. Note that we do not include radiation and baryon in this investigation. 
However, if we add all constituents in our investigation, there will be the unstable fixed 
points in which the added constituents dominate. Properties of these fixed points are similar 
to those of fixed point (a). The eigenvalues of these fixed points can be generalized as 
[ —3, 3(1 + Wj)/2, 3(1 + oji) ], where Wj is the equation of state parameter of the dominated 
constituent i. 

The fixed points (b) correspond to the three-form field dominated solution and give 
w tot = — 1- Therefore, three- form field behaves as dark energy. They give the eigenvalues 
(—3,0,-3/2). Since one of the eigenvalues is zero, we can say nothing about the stability 
of the fixed points from the linear analysis. We have to use the second order perturbations 
and specify the potential to analyze the stability of these fixed points. However, we will 
use the numerical method to analyze the property of the fixed points. In the numerical 
method, we check whether the small deviations of the dynamical variables are convergent 
and also check whether their evolution provides the standard evolution of the universe. Note 
that the standard evolution of the universe we use in this paper means that there are the 
radiation- matter- and dark energy-dominated periods respectively. For the exponential 
potential, V = Vbexp(— rjx), we find that the fixed point (bl) is stable for 77 > and the 
fixed point (b2) is stable for rj < 0. For inverse power law potential, V = Vox~' n and 
Gaussian potential V = Vq exp(— rjx 2 ) both of the fixed points (b) are stable for rj > — 1 and 
rj > respectively. These properties of the fixed points (b) are consistent to our analytic 
investigation that the potential should be runaway and it is compatible with the second order 
perturbation analysis in [20]. 

The fixed points (c) provides the cosmological constant behavior since it corresponds 
to A = for runaway potential. This is the special case of the runaway potential since it 
yields the constant energy density of dark energy due to V = Vq. We note that it may 
not be completely correct to mention that it is runaway potential when it is constant. This 
constant potential leads to the stable fixed point with y 2 + z 2 = Qx = 1> ^tot = — 1 an d 
wx = constant = — 1. This fixed point corresponds to the extremum fixed points in [20]. We 
avoid the call this fixed point as extremum fixed point since it is constant potential. We also 
note that there is the extremum fixed point for Gaussian potential with x = 0,y = 0, z = 1 
and A / 0. This fixed point will be stable if rj < 0. However, it leads to the dark energy 
dominated period in the past. 

Cosmological constant model can be promoted as a subclass of three-form dark energy 
model. This is a useful feature of dark energy from three-form field. However, it cannot 
solve the coincidence problem which state that "why dark energy is comparable to matter 
nowadays" . This is due to the fact that there are no the stable fixed points at which Qx /^c — 
7/3 where Qx and f2 c are the energy density parameter for dark energy and dark matter 
respectively. Next section, we will investigate an probability for solving this problem by 
introducing the interactions between dark energy and dark matter. 
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3 Coupling Three-form Field with Dark Matter 



According to observations, today the dark matter density is close in value to the dark energy 
density. This leads to the so-called coincidence problem because their evolutions are consid- 
erably different throughout the universe history. Is it coincidence that their energy densities 
are of the same order? There have been many models proposed to explain the coincidence 
problem and it has been found that an uncoupled three-form field cannot solve the coinci- 
dence problem. If dark matter is capable of decaying into dark energy, the explanation of 
the similarity of their energy densities may be made. This introduces the coupling between 
three-form field dark energy and dark matter. In order to alleviate the coincidence problem 
we expect this coupling to lead to an accelerated scaling attractor solution 

^ darkenergy = 0(1) and a > 0. (3.1) 

^ Mark matter 

The existence of the coupling can be represented by the modified continuity equations 

p c = -m Pc - Q, (3.2) 
px = -3H( Px +px) + Q, (3.3) 

where p c stands for energy density of cold dark matter and Q is the energy transfer between 
dark energy and dark matter 

Q > =4> dark matter — > dark energy 
Q < => dark energy — > dark matter, 

while the background baryons and radiation still satisfy 

p b = -3Hp b , (3.4) 

p r = -3H(p r +p r ), (3.5) 

where b stands for baryon and r radiation. The explicit form of the Einstein equations (2.3) 
and (2.4) including baryons and radiation become 



H 2 = - 



K 



^(X + 3HX) 2 + V(X) + Pb + p r + Pc 



(3.6) 



K 2 



H = — — (V,x X + pf, + p r + p r + p c ). (3.7) 

From equation (3.3), using (2.15) and (2.16) we obtain 

X + 3HX + 3HX + V, x = -■ — • (3.8) 

X + 3HX 

The interaction term we added in equations (3.2) and (3.3) will modify the dimensionless 
variable equations (2.19) and (2.20) as 

y' = 7 - \\{x)z 2 (xy - ^|) + ^(w 2 + V + v 2 )y, (3.9) 

v > = -11 _ l v u + A(x)xz 2 - w 2 - -u 2 - v\ (3.10) 
v 2 3 



- 8 - 



where 

7= F ■ KQ . (3.11) 

>/6(Jf + 3HX)H 2 

Coupling forms in which we use to consider in this paper are phenomenological. The coupling 
forms are introduced by taking into account the fact that it should be expressed in the 
covariant form due to the ability to compute the cosmological perturbations. According to 
the quintessence model, there are three simple coupling forms. The first coupling is Q oc (f>p c , 
where <j) is the quintessence field. This coupling form is motivated from the scalar-tensor 
theory [10]. Unfortunately, by considering the general form of Lagrangian of scalar field 
with this coupling form, there are no fixed points giving the matter- or radiation-dominated 
periods if dark energy dominates at the present time [21]. Analogous to this coupling form of 
the scalar field and its simple kinetic term, we introduce the coupling form of the three-form 
field by replacing <j) with X + 3HX. Note that there is no theoretical motivation like the 
quintessence models. 

The second class of the coupling is Q oc H{a ( pp ( j > + a c p c ), where p^ is energy density of 
the scalar field and and a c are dimensionless constants. The motivation of this coupling 
form does not come from the fundamental theory but mathematical calculation in which the 
energy density ratio of dark matter to dark energy is constant at a stable fixed point [30]. 
In this coupling form, there is no the problem like the previous one but it gives a negative 
energy density of dark energy at matter- and radiation-dominated period [22]. By considering 
au = case with exponential potential, the coincidence problem cannot be solved since there 
are no stable fixed points Q^/Qc ~ 0(1) [12]. Note that the transfer rate of this coupling 
form is determined by the expansion rate of the universe rather than by the local quantities 
associated with particle interactions. Therefore, it suggests that it is not easy to describe 
this coupling form in terms of physical interaction of particles. 

The third coupling form is introduced in order to obtain the physical motivation asso- 
ciated with particle interactions. The similar form of this coupling has been use in reheat- 
ing, decay of curvaton and dark matter into radiation [12]. This coupling takes the form 
Q oc (r^/9^ + r c p c ), where is a constant. This coupling form does not depend on the ex- 
pansion rate of the universe. T acts as a local transfer rate which is similar to the decay rate 
of the particle interaction in collision approach. This form is quite complicated because one 
cannot eliminate the Hubble parameter from the dynamical equations of the dimensionless 
variables. Similar to second coupling form, it encounters the negative energy density in full 
parameters consideration (r^, 7^ and T c 7^ 0) while it cannot solve the coincidence problem 
in the case of = [23] . 

It is worthwhile to investigate whether interacting three-form field can yield the solution 
of the coincidence problem. In order to study the three-form interaction with dark matter 
phenomenologically, we introduce the covariant forms of the coupling for three-form field as: 



(I) = 7273k/3T; (c) (l/24a 3 ) e^F upa ^, 

(II) Q" = aHTT (c) u^ 



r(c)' 

(in) Q^ = rr; (c) < 

where a, /3 and T are constant, is a 4- velocity, TL c ^ = —p c is the trace of the energy- 
momentum tensor for cold dark matter. is the covariant form of energy-momentum 
transfer corresponding to the covariant form of the conservation equation [31] 

V./Zf = Qf, (3.12) 
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with 

X)Q? = 0, (3.13) 

i 

where i is a component of the universe. The zero component of the covariant conservation 
equation (3.12) corresponds to the conservation equations of the energy density (3.2) and 
(3.3). This turns the relation between covariant 4- vector and the energy transfer Q as 

Q = Qo- (3.14) 

It is important to note that the coupling model (II) is not the expression of the covariant 
form due to the appearance of H. However, the covariant forms are introduced in order to 
provide the ability to compute the perturbations [12] and the coupling model (II) can provide 
the ability to compute the perturbations which has been performed in [32]. Therefore, we 
allow the appearance of H to arise in the coupling model (II). 

Next we will consider the properties for each coupling form in detail by finding the 
fixed points and then investigating their stability. For simplicity, we will use the exponential 
potential (V = Voe - ^), the Gaussian potential (V = Voe~ vx ), and the inverse power law 
potential (V = Vqx~ v ) in our investigation. The evolution behavior is also investigated in 
order to compare with the standard evolution of the universe. Conveniently, we consider only 
the late time fixed points since they are relevant to the dynamics of dark energy. We also 
ignore the effect of baryon contribution in qualitative analysis because it has a tiny effect on 
the dynamical evolution comparable to other constituents. However, in the computational 
simulation of dynamical evolution, we add up all constituents in our consideration. 

3.1 Coupling model (I): = ^/2/3k^TL c) (l/24a 3 ) e v P^F upa ^ 

For this covariant form, the interaction term that satisfies equations (3.2) and (3.3) can be 
expressed as Q = Qo = y / 2/3K(3p c (X + 3HX) and the interaction variable becomes 7 = (3v 2 . 
The autonomous system in (3.9) and (3.10) can be rewritten as 

y > = _ *-\{x)z\xy - y|) + \{w 2 + \u 2 + v 2 )y, (3.15) 
v ' = -f3 yv - ^-v[l + X(x)xz 2 - w 2 - \u 2 - v 2 ]. (3.16) 

The crucial properties of the fixed points for this interaction are summarized in Table 2. 
Next we will discuss the results of each fixed point in detail. 

• Fixed points Al and A2 

These fixed points correspond to z = 0. Since o; t ot depends on parameter /3, it is easy to show 
that accelerating phase occurs when > \/3/2. Notice that the phantom phase w to t < — 1 
is impossible because this will cause t>* to be imaginary. The range of (3 is constrained by 
the existence condition |/3| < 3/2. The stability of these fixed points can be analyzed by 
considering the eigenvalues from the first order perturbation around each fixed point of the 
autonomous system (equation (2.18), (3.15), and (3.16)). These fixed points will be stable 
if the real part of all eigenvalues are negative. For these fixed points, the eigenvalues do not 
depend on the potential and read (—3, (9 — 4/3 2 ) /6, (9 — 4/3 2 ) /3). This yields the stability 
condition |/3| > 3/2. Since the stability condition is not compatible with the existence 
condition, we conclude that Al and A2 are not stable. We are not interested in these fixed 
points. 
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Table 2. The properties of the fixed points for the coupling model (I) with the exponential potential. 
Note that x = y/2/3y at the fixed points. 



• Fixed points Bl and B2 




Figure 1. This shows the evolution of the dynamical variables with different interaction parameter 
(/?) at fixed point B, = +1. The red line represents the energy parameter of dark matter (£l c ), 
the blue line represents the energy parameter of dark energy (fix) and the green line represents the 
energy parameter of radiation (fl r ). In this plot, we use the exponential potential with A = 1 and 
P = 0.5 for (a), p = for (b), and /3 = -0.5 for (c). 

These fixed points correspond to dark energy dominated point and require z = 0. We cannot 
use these fixed points to solve coincidence problem because Qx/^c 1S n °t order of unity. 
However, it is instructive to investigate the effect of the coupling term on the dynamical 
evolution. The most important thing of this coupling form is that one can reproduce the 
non-interacting approach. There are no boundaries for the existence condition. For the 
stability condition, the eigenvalues do not depend on the potential and can be expressed as 
(—3, 0, —3/2-/3). Thus, it guarantees that the fixed points are not stable for j3 < —3/2. In the 
case of P > —3/2, due to the zero eigenvalue, we have to use the second order perturbations 
for analyzing the stability condition as seen in [20] or use numerical method as we use here. 

In numerical method, we check whether the evolution of the dynamical variables is 
convergent and provides the consistent behavior. For the exponential potential, the result 
is that if y* = 1, the potential parameter needs to be positive (rj > 0) and if y* = —1, the 
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potential parameter needs to be negative (77 < 0). In other words, it is only allowed for 
7?y* > 0. From the first term in equation (3.16), the transfer energy also depends on the 
sign of y*. The result is that /3y* > corresponds to transfer of dark matter to dark energy 
and f3y* < corresponds to transfer of dark energy to dark matter. However, if we carefully 
consider the evolution behavior of the dynamical variables, we will see that the interaction 
parameter (/3) cannot take more value. This is due to the fact that the more value of the 
interaction parameter, the more change of the behavior of the variables at the early time. 
Figure 1 shows the effect of the interaction on energy density parameter evolution. We use 
the exponential potential with /3 = 0.5 in Figure 1(a), /3 = in Figure 1(b), and /3 = —0.5 
in Figure 1(c). For these fixed points, they are allowed for both energy transfer from dark 
energy to dark matter and dark mater to dark energy. 

The evolution for the Gaussian potential is similar to the exponential potential. The 
difference is that it is only stable for the rj > case. This result agrees with the analytic 
version in [20]. For the inverse power law potential, the evolution is stable if 77 > — 1. It is 
consistent to the qualitative analysis as we have analyzed before. 




Figure 2. (a): This figure shows the region of stability (red, shaded) and existence (blue, shaded) for 
the exponential potential in the (/3, rj) parameter space. The violet region indicates the compatible 
region of two conditions. Dashed line represents a point at which the energy parameter ratio of dark 
energy to dark matter is 7:3. (b): This figure shows the evolution of the dynamical variables for 
the exponential potential. We use /3 = —2.13 for this simulation. The red line represents the energy 
parameter of dark matter (f2 c ), the blue line represents the energy parameter of dark energy (fix) 
and the green line represents the energy parameter of radiation (fi r ). 

• Fixed points CI and C2 

These fixed points correspond to dark energy dominated period at which u; tot = — 1. 
Moreover, these fixed points provide us the scaling regime. However, by considering the 
allowed region of the parameters for the stability and existence conditions, it cannot provide 
the well-behaved evolution as we will show now. 

By using the same step, the plots of stability and existence region for the exponential 
potential are shown in the Figure 2(a). The dashed line, which shows the points at which 
the ratio of dark energy to dark matter is 7:3, is in the region that satisfies two conditions. 
In order to compare the result of this model with the observational data and solve the 
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coincidence problem, f3 and r\ must satisfy the relation 



2\/6/3 3 
45 - 14/3 2 



(3.17) 



This relation suggests that 77 is divergent at ft = ±^45/14 as shown in Figure 2(a). This 
divergence also exists in the Gaussian potential and the inverse power law potential. From 
Figure 2(a), it is found that we will solve the coincidence problem by using these fixed points 
as long as |/3| > y 45/14. Due to the more value of the coupling parameter, it is obvious 
that the evolution at the early time will not satisfy the standard evolution of the universe as 
shown in Figure 2(b). We will see that there is no matter-dominated period in the evolution. 
The disappearance of matter-dominated period is also found for the Gaussian potential and 
the inverse power law potential. 

Finally, it is important to note that there is the cosmological constant fixed point, A = 0, 
like non-interacting case. The interaction terms do not affect this fixed point because the 
interaction terms depend on v and v = at this fixed point. Properties of this fixed point are 
similar to the non-interacting case. The difference is only that the matter does not completely 
dominate at matter-dominated period. 



3.2 Coupling model (II): 



aHT T , .vf 

t(c) 



The interaction term that satisfies equations (3.2) and (3.3) of this coupling model can be 
expressed as Q = Qo = aHp c and the interaction variable becomes 7 = av 2 /2y. The 
autonomous system in (3.9) and (3.10) can be rewritten as 



— av v\l + X(x)xz 2 

2 2 1 { ' 



3) + 2^ +3" + v )2/> 



w 



-u 
3 



(3.18) 
(3.19) 



The properties of each fixed point are summarized in Table 3. Note that there are the other 
fixed points at which the interaction parameter is a = —3. It is not an interesting fixed point 
due to the much more value of a. As we discussed in the previous interaction, the more value 
of a, the more deviation from standard evolution of the universe in the past. Moreover, this 
value provides the non-dynamical evolution of the dark matter energy density, p c = 0. This 
behavior is not accepted for the observational data. 
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Table 3. The properties of the fixed points for the coupling model (II) with the exponential potential. 
Note that = -\/2/3y* at the fixed points. 
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• Fixed points Al and A2 

These fixed points correspond to z = and yield the phantom phase when a < — 3. 
However, this phase does not exist since v is not real number. The existence condition 
for these fixed points is the same for all potentials. This condition can be expressed as 
—3 < a < 0. To find the stability condition, we do the linear perturbations around the fixed 
points. The eigenvalues of these fixed points read (—3, a + 3, a + 3). We will see that they 
do not depend on the potential form. Thus the stability condition is expressed as a < — 3. 
This condition conflicts with the existence condition for all potentials. This means that these 
fixed points are not stable for all potentials. We are not interested in these fixed points. 

• Fixed points Bl and B2 




(a) (b) 

Figure 3. (a) shows the evolution of the dynamical variables for y„ = +1 fixed point with the 
exponential potential. We use a = —0.2, rj = 1.0. The red line represents the energy parameter of 
dark matter (f2 c ), the blue line represents the energy parameter of dark energy (Ox) and the green 
line represents the energy parameter of radiation (f2 r ). (b) shows the evolution behavior of the total 
equation of state parameter, w to t, by using the same parameters with (a). 

These fixed points correspond to dark energy dominated point without the effect of 
potential, z = 0. These fixed points cannot be the points at which the coincidence problem is 
solved because of Qx = 1- However, we will investigate the properties of these fixed points. 
These fixed points satisfy the existence condition for all points in parameter space (ct,r]). 
This condition is valid for all potentials. 

Considering the stability condition, the eigenvalues of the linear perturbation system 
do not depend on the potential form and read (—3,0, —3/2 — a/2). As we discussed in the 
coupling model (I), a zero of the eigenvalues indicates that one cannot find the stability 
condition of these fixed points at the linear order perturbations. We check this by using a 
numerical method as we have done in the coupling model (I). Interestingly, it has similar 
properties to the coupling model (I). From the above eigenvalues, the stability condition is 
a > —3. The crucial property of this coupling depends on the sign of a. For the case of 
a < with the exponential potential, = +1 fixed point is stable when rj > 0; conversely, 
y* = — 1 fixed point is stable when rj < 0. In other words, they will be stable if r/y* > 0. For 
the Gaussian potential and the inverse power law potential, both y* = ±1 fixed points will 
be stable if r\ > and rj > — 1 respectively. 
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Figure 4. (a) shows the evolution of the dynamical variables for = +1 fixed point with the inverse 
power law potential. We use a = —0.2, rj = 1.0. The red line represents the energy parameter of 
dark matter (fi c ), the blue line represents the energy parameter of dark energy (Qx) and the green 
line represents the energy parameter of radiation (f2 r ). (b) shows the evolution behavior of the total 
equation of state parameter, uj to t, by using the same parameters with (a). 



For the case of a > 0, corresponding to energy transfer from dark matter to dark 
energy, with the exponential potential, we find that the evolution from numerical simulation 
encounters a singularity at the matter-dominated era. This singularity can be seen from 
the interaction term of dynamical equation (3.18) when !) / and y = 0. However, this 
singularity does not exist in the case of a < as shown in Figure 3. From this figure, we will 
see that £lx does not completely vanish at the matter-dominated period. This behavior also 
exists in this class of coupling form for the quintessence model as seen in [33] and references 
therein. The physical interpreting of this behavior is that the energy density of dark matter 
is higher in value in the past than the present value of the energy density. This leads to the 
point that v — >■ 1 and y — > 0. Moreover, this coupling form also provides a clue to violate 
the second law of thermodynamics similar to the quintessence model [34]. Note that this 
singularity also exists in the Gaussian potential and the inverse power law potential when 
a > 0. 

From Figure 3(b), we will see that uj to t does not vanish at the matter-dominated period, 
N ~ —3, since O c does not completely dominate. It is also found that u) tot oscillates at the 
radiation-dominated period. This effect also occurs in all coupling models we investigate. To 
see how this effect can occur, we will return to the general form of the total equation of state 
parameter, 

wtot = -1 + {irU 2 + 7 c w 2 + ^bW 2 ~ \xz 2 ). (3.20) 

From this general form, we will see that wtot = at the radiation-dominated period if z 
completely vanishes. However, according to simulation computing, z does not completely 
vanish. It oscillates in order of 10~ 9 . As we have mentioned before, a; is a large value at early 
time. Thus, the large value of x will factor the last term in equation (3.20) to oscillate in 
order of the huge value. However, for the inverse power law potential, x will be eliminated by 
A = 7]/x and the oscillating behavior will disappear as shown in Figure 4(b). From this figure, 
it is found that the oscillating behavior disappears and wtot ~~ ^ w r at radiation-dominated 
period as inferred from equation (3.20). This is the advantage of inverse power law potential. 
Thus, the inverse power law potential is favored by our numerical calculation and it implies 
that it is the most useful potential among all potentials we have investigated. This argument 
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will hold only in the background analysis and may not be valid when the perturbation analysis 
is taken into account as suggested in [20]. We leave this work in further investigation since 
it is beyond the purpose of our study. 

The cosmological constant fixed point, A = 0, also exists in this coupling model. The 
oscillating behavior also disappears since the last term in (3.20) completely vanishes. The 
well-behaved evolution of the dynamical variables is allowed only for a < 0, corresponding 
to the energy transfer from dark energy to dark matter. The physical interpreting of this 
property is the same as we have mentioned above. 

3.3 Coupling model (III): Q» = VT^u^ 1 
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Table 4. The properties of the fixed points for the coupling model (III) with the exponential potential. 
Note that cc* = ^2/3y* at the fixed points. 



The interaction term that satisfies equations (3.2) and (3.3) of this coupling model can 
be expressed as Q = Qq = Tp c . We cannot eliminate the Hubble parameter from the equation 
system of the dimensionless variables. We have to introduce new dimensionless variable for 
this coupling model. This variable has been introduced as s = Hq/(H + Hq) [12]. This form 
of this variable satisfies the compactness property but provides the singularity in the equation 
system at s = 1. In our consideration, we ignore the compactness property and choose the 
form of this variable for the simplicity of calculation. This new variable is expressed as 

s = § , (3.21) 

where Hq is the Hubble parameter at the present time. Due to the fact that H oc 1/t, it 
implies that s — > at the early time and s = 1 at the present time. The interaction variable 
becomes 7 = jsv 2 /2y, where 7 = Y /Hq. The autonomous system in (3.9) and (3.10) can be 
rewritten as 

/ jsv 2 3 2 2 2 2i 1 v6,/ \ 2 /, no i 

y = -7, ~V A(sc)xz - -f b w - 7 r u -j c v \ + -rrK x ) z i I 3 - 22 ) 

2y 2 2 

v' = -l^L - ^v[j c + \(x)xz 2 - j b w 2 - 7 r n 2 - 7 C ?; 2 ], (3.23) 
and the equation of new dimensionless variable can be expressed as 

s' = ~s[\{x)xz 2 - lb w 2 - lr u 2 - ~f c v 2 }. (3.24) 
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As we have seen, s does not need to be stable at the present or at the early time due to the 
fact that it scales as 1/H. It is well known that the Hubble parameter is not stable in all era 
of the late-time history of the universe (except inflationary era). However, we will consider 
in both cases, considering a fixed point of s and ignoring it. The properties of each point are 
summarized in Table 4. 

• Fixed points Al, A2 

These fixed points are (x,y,v,w,u, s) = (0, ^—5*7/3, ±yl + s*7/3, 0, 0, s*). This is the 
fixed points at which s does not need to be a fixed point. If we require the fixed point of s, it 
turns out that s = 0. These fixed points can be reduced to (x, y, v, w, u, s) = (0, 0, ±1,0, 0, 0) 
and corresponds to dark matter dominated solution. Its eigenvalues are (—3,3,1/2) for all 
the potentials. Therefore, these fixed points are not stable. However, if we return to ignore 
the fixed point of s, the eigenvalues will be (—3, 3 + 7s, 3 + 7s). Thus, the stability condition 
can be written as 7s < —3. This condition conflicts with existence condition, —3 ^ 7s ^ 0, 
and implies that these fixed points are not stable. 




(a) (b) 

Figure 5. (a) shows the evolution of the dynamical variables for t/„ = +1 fixed point with the inverse 
power law potential. We use 7 = —0.2, r\ — 1.0. The red line represents the energy parameter of 
dark matter (f2 c ), the blue line represents the energy parameter of dark energy (Ox) and the green 
line represents the energy parameter of radiation (f2 r ). (b) shows the evolution behavior of the total 
equation of state parameter, cj to t, by using the same parameters with (a). 



• Fixed point Bl, B2 

These fixed points are (x,y,v,w,u, s) = (±y/2/3, ±1, 0, 0, 0, s*). They correspond to three- 
form dominated solution and their properties are the same for all potentials. It gives the total 
equation of state parameter like the cosmological constant, uj tot = —1. The eigenvalues of 
these fixed points are (—3, 0, (— 3— 7s*)/2). Therefore, we can infer nothing about its stability 
from the linear analysis because of the eigenvalue 0. However, the last eigenvalue suggests 
us that the stability condition requires at least 7s* > —3. If this condition is satisfied, the 
stability condition will be investigated by using numerical method as done in the previous 
coupling models. It turns out that these fixed points are stable. The crucial properties of 
these fixed points are similar to the coupling model (II). The well behavior of the energy 
density evolution is allowed only for negative sign of interaction parameter, 7 < 0. In other 
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words, it is allowed for energy transfer from dark energy to dark matter only. The physical 
interpretation of this phenomenon is the same as mentioned in the coupling model (II) . 

The advantage of this coupling form is that it has a tiny effect on the evolution behavior 
at the matter- and radiation-dominated periods as shown in Figure 5. From this figure, the 
evolution behavior at the matter- and radiation-dominated periods is close to the standard 
evolution of the universe. This feature is explained by considering the interaction terms in 
equations (3.22) and (3.23). These interaction terms effectively vanish since s approaches 
zero during matter- and radiation-dominated periods. They also have a tiny effect at the 
future time since v — > for these fixed points. This keeps the stability of these fixed points 
for all variables except s even though s still evolves. 

• Fixed point CI, C2 

These fixed points are (x, y, v, w, u,s) = ( 0,0, s*), where x* = y/2/3y*, < y* < 1, 



The most important properties of these fixed points are related to the second relation (3.26). 
This relation yields T = —3H. It provides not only the more value of the interaction param- 
eter but also the non-dynamics of energy density of dark matter due to equation (3.2). This 
gives the similar results to the case of a = —3 in the coupling model (II) by replacing 7s* 
with a. The difference is that 7s* varies with time and approaches zero at the early time 
while a is constant. However, the evolution behavior is not significantly different since the 
matter will run out before it dominates in the past due to the more value of 7s* at the fixed 
points. 

For the case of s to be a fixed point, it is found that s* = and this yields the infinity of 
the coupling parameter, 7 = — 3/s* — > 00. If we take a small perturbation by setting = e 
where e is a small value, it turns out that this fixed points are not stable and evolve to fixed 
points B. 

Finally, we note that there is the cosmological constant fixed point, A = 0, similar to 
coupling model (I) and (II). The difference is that the evolution behavior are tinily affected 
from the interaction since s — > at the matter-dominated period. We also found that one 
can set s as any value at this stable fixed point due to v = fixed point. 

4 Extension to Another Coupling Model 

We will begin this section by summarizing the results of three coupling models. 
• Coupling model (I) 

There are three classes of all fixed points. The Al and A2 fixed points are not stable. For the 
Bl and B2 fixed points, they are stable but one cannot use them to solve coincidence problem 
since they are the fixed points at which y = ±1. However, it is instructive to consider this 
fixed points in order to study the effect of the interaction because it can be reduced to the 
non-interacting approach. The energy can transfer from dark energy to dark matter and 




(3.25) 



and 



7 s * 



= -3. 



(3.26) 
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dark matter to dark energy. Finally, for fixed points CI and C2, the coupling parameter is 
allowed only in region with large value. Therefore, there is no matter-dominated period for 
these fixed points. 

• Coupling model (II) 

Actually, there are three classes for this coupling model, but we show only two classes because 
the missing one provides the much more value of the coupling parameter, a = —3. For the 
fixed points Al and A2, they are not stable like the fixed points Al and A2 in the coupling 
(I). The coincidence problem cannot be solved by using the fixed points Bl and B2 because 
of y = dbl. The energy transfer is allowed only from dark energy to dark matter. According 
to the evolution behavior, we found that the inverse power law potential is suitable for the 
interacting three-form field. 

• Coupling model (III) 

The properties are very similar to the coupling model (II) although there are more variables 
in phase space. They are similar when we take 7s* — > a. This means that all properties are 
the same as the coupling model (II) except that the effective coupling 7s* now varies with 
time. 7s* is equal to —3 at the fixed points but not necessary hold for other points in phase 
space. This provides the vanishing value of effective coupling at the matter-dominated period 
because of s — >■ 0. However, the evolution behavior of these fixed points is not compatible 
with observation since there is no matter-dominated period. In other words, the matter will 
run out before it dominates in the past due to the more value of 7s* at the fixed points. 

The cosmological constant fixed point, A = 0, also exists in all coupling models we have 
investigated. It is the stable fixed point with px = constant and wj = —1. Therefore, the 
coincidence problem will not be solved by this type of the fixed point. 

From summary above, there are no coupling models which provide the possibility to 
solve the coincidence problem. We now introduce the another coupling model in order to 
solve it. This coupling model can be written in the covariant form as 

where 7 is the coupling parameter for this model. For this covariant form, the interaction 
term that satisfies equations (3.2) and (3.3) can be expressed as Q = Qo = V6T(X+3HX)/k. 
One cannot eliminate H from the dimensionless equations as found in coupling model (III). 
We introduce new dimensionless variable s = Hq/H as done in coupling model (III). We 
note that it is possible to eliminate H from the dimensionless equations by introducing the 
covariant interaction form as 

1 

«" = -^-^ mF -"-""- < 4 - 2) 

However, it encounters the singularity at the radiation-dominated period. We will discuss 
this issue later. The interaction term in equation (4.1) yields 7 = 7s 2 , where V = Hq^. The 
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dynamical equations associated with the interaction term become 



y' = is 2 - ^X(x)z 2 {xy - J |) + ^(w 2 + ^u 2 + v 2 )y, (4.3) 

v , = _j^y _ 3 + 2 _ w2 _ 4^ 2 _ ^ 

v 2 3 

3 

s' = --s[X(x)xz 2 - j b w 2 - JrU 2 - jcV 2 }- (4.5) 

Note that the dynamical equations for the interaction (4.2) are obtained by ignoring the 
equation (4.5) and setting the interaction parameter as 7s 2 — > 7. From the first term in 
equation (4.4), we will see that the energy density of dark matter cannot vanish. This 
guarantees that there are no stable fixed points at which J7 C = at the present time. It is 
also found that this term will provide the singularity at radiation-dominated period since 
v — > for the interaction (4.2). However, for the interaction (4.1), the singularity can be 
avoided since s — > at the radiation-dominated period. This mechanism is similar to a 
mechanism in the coupling model (III) in which the effect of the interaction is very small 
at early time. Therefore, the interaction in (4.2) is not interesting for us and we will not 
consider it. The properties of all fixed points are summarized in Table 5. Conveniently, we 
introduce the new parameters as 

1/3 



A = 97^ + 3^-3 + 9(752)2 , (4.6) 



B^-( 7 sl) 2 -3V6^. (4.7) 



• Fixed points A1-A6 

These fixed points correspond to = 0. The constraint equation for these fixed points can 
be written as 

yl -y*- \isl = 0. (4.8) 

Three solutions of for this equation are expressed in Table 5. The conditions for stability 
and existence of all A1-A6 do not depend on potential form. These two conditions are not 
compatible each other. We are not interested in all these fixed points. 

• Fixed points B1-B4 

These fixed points correspond to the solutions with ^ 0. One of the conditions for these 
fixed points is 

Ax*;? 2 - 7 b w 2 - 7 r n 2 - j c v 2 = 0. (4.9) 

This condition leads to co t ot = = — 1 and also provides the fixed point s* 7^ 0. It is 
possible to set s* = 1 as the fixed point at the present time to get H = Hq. The condition 
(4.9) can be rewritten as 

yl - l£y* - f 1 + yJl^Y J = and vl = -psly*. (4.10) 
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Table 5. The properties of the fixed points for the coupling model (IV). Note that x* = y / 2/3y* at 
the fixed points. We set s* = 1 for convenience. To restore s* into the table, one just replaces 7 by 
7s*. For fixed point B, solutions of y* depend on the potential form. In this table, we show only for 
the solution of exponential potential. For the other potentials, we will discuss in the description of 
each potential. 



The solution of is strongly dependent on the potential form. There are two solutions 
for the exponential potential and inverse power law potential and three solutions for the 
Gaussian potential. From this condition, we will see that it is possible to get the fixed points 
at which vl/yl ~ 0(1) with = 1. This allows us to solve the coincidence problem since 
Oo/Oyo = ^*/(y* + z 1)- Our investigation of these fixed points will be separated into three 
parts depending on the potential. 

• Inverse power law potential 

For inverse power law potential, A = r)/x* = v/3/2r//i/*. By substituting A into equation 
(4.10), the solutions of this equation are 

2 7 s 2 + 2 7 s 2 ?? ± V36i? 2 + (-2 7 s 2 - 275^)2 

y*± = 7. • (4.11) 

077 

The conditions of existence and stability are shown in Figure 6. The minus solution 
corresponds to the Figure 6(a), and plus solution corresponds to the Figure 6(b). By using 
the solutions in equation (4.11), one can find the relation of 7, 77, s* and Vl c as 

„ = McM) 2 (412) 

V 9^2 -4(1 -a,)^- ^ • > 
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Figure 6. This figure shows the region of stability (red, shaded) and existence (blue, shaded) in the 
(7, 77) parameter space for the inverse power law potential V = Vox~ v . The violet region indicates the 
compatible region of two conditions. Dashed line represents a point at which the energy parameter 
ratio of dark energy to dark matter is 7:3. Figure (a) and (b) correspond to y*_ and y*+ solution in 
equation (4.11). Note that we set s* = 1 in this plot. However, we can restore s* into this plot by 
setting 7 — > 7 si. 



The characteristic properties of this relation are found in the Figure 6. From this 
relation, the singularity of the curve £lx /^c m Figure 6 can be written as 



For f2 c = 0.3, we get the singularity at 7 s 2 ~ ±0.538 as shown in Figure 6. From Figure 
6 and relation (4.12), It is also found that the stability condition can be written as 

i*i > A: (4 ' 14) 

This condition also suggests us that the strength of the interaction depends on the value of 
Q c . If fi c increases, the lower bound of this condition also increases. Thus £l c cannot take 
large value since interaction will be too large. 

For the evolution behavior, we set 7 = —0.5, s* = 1 and choose Q c ~ 0.25. This 
leads to r\ = 167 2 /3(— 3 + I67 2 ) « 1.333. By using this parameter set, the evolution of 
energy density parameters and Wtot are shown in Figure 7. It is found that one can solve the 
coincidence problem by using this coupling model because there are the stable fixed points 
at which f^x/^c ~ 7/3 and there are the matter- and radiation-dominated periods with the 
same parameters. Note that, in the simulation, we have to set a small deviation from the 
fixed points in order to get the dynamics of the variables. Therefore, the stable fixed points 
will be slightly deviated from the fixed point we set as found in Figure 7. 

• Exponential potential 
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Figure 7. (a) shows the evolution of the dynamical variables in extension the coupling model with 
the inverse power law potential. We use 7 = —0.5, rj = 3 (_3 ( | 7 167 a-) and s* = 1 for this simulation. 
The red line represents the energy parameter of dark matter (J7 C ), the blue line represents the energy 
parameter of dark energy (fix ) and the green line represents the energy parameter of radiation (fi r ). 
(b) shows the evolution behavior of the total equation of state parameter by using the same parameters 
with (a). 
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7 7 
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Figure 8. This figure shows the region of stability (red, shaded) and existence (blue, shaded) in 
the (7,7?) parameter space for the exponential potential V = Voe~ vx . The violet region indicates the 
compatible region of two conditions. Dashed line represents a point at which the energy parameter 
ratio of dark energy to dark matter is 7:3. Figure (a) and (b) correspond to y+ and y_ solution 
respectively. Note that we set s* = 1 in this plot. However, we can restore into this plot by setting 
7 -> is 2 * ■ 



For the exponential potential, A = 77 and the solutions of corresponding to the equation 
(4.10) are 



y*± = 5 • (4.15) 
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The plus solution corresponds to the fixed points Bl and B2, and minus solution corresponds 
to the fixed points B3 and B4. We show the existence and stability conditions by using 
numerical method as seen in Figure 8. From this figure, only y* + solution is consistent while 
there are no existence regions for solution. Thus, the interesting fixed points are Bl and 
B2 which have the same properties. For S]x/U c = 7/3, the stable condition allows only the 
region in which | | ^ ^81/280 « 0.538. 

The calculation and analysis are performed in the same step as done in the inverse 
power law potential. The results are similar to the inverse power law case. The difference is 
that the evolution behavior for exponential potential will encounter the large oscillating of 
cJtot a t radiation-dominated period. We will not show this evolution behavior here. 




-2 -1 1 2 -4 -2 2 4 



t y y 

(a) (b) (c) 

Figure 9. This figure shows the region of stability (red, shaded) and existence (blue, shaded) in 
the (j,t]) parameter space for the exponential potential V = Voe~ vx . The violet region indicates the 
compatible region of two conditions. Dashed line represents a point at which the energy parameter 
ratio of dark energy to dark matter is 7:3. Figure (a), (b), and (c) correspond to three solutions of 
equation (4.10). 



• Gaussian potential 
For the Gaussian potential, A = 2r/x* = y^8/3r]y^. Equation (4.10) become 

y* ^ i - 7) — = °- 4 - 16 ) 

3 2ryy* 

There are three solutions of this equation. We will not show the exact solutions here 
for convenience because they are lengthy. The region plot of these three solutions are shown 
in Figure 9. The evolution behavior is similar to the inverse power law and exponential 
potential case. Thus, we will not show this simulation here. 

Finally, it is important to note that there is no cosmological constant fixed point, A = 0, 
similar to three coupling models we have investigated. This can be seen from the interaction 
terms of this coupling model, (4.3) and (4.4). These terms do not vanish when v — > 0. 
Therefore, it yields the results equivalent to the fixed points A which = and vl = 
— 27s^/3y*. It is not a stable fixed point for this coupling model. 
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5 Conclusions and Discussions 



The possibility of the three- form field to be dark energy is investigated. We begin this work 
by summarizing the results of three-form field dark energy. By analyzing the dynamical 
equations of the system, we have found that the potential for three-form field should be 
runaway. We have also found that three-form field with proper potential can act as dark 
energy. However, non-interacting three-form field cannot provide the possibility to solve the 
coincidence problem since there are no stable fixed points at which Qx/^c ~ 0(1). Analogous 
to scalar field dark energy, we examine the interacting three-form field with dark matter and 
find the possibility to solve the coincidence problem. We introduce three covariant coupling 
forms which imitate idea from scalar coupling forms. Although the coupling model (II) can 
not be expressed in the covariant form due to the appearance of H, we can compute the 
perturbations for this coupling model as seen in [32]. Therefore, we allow the appearance 
of H to arise in this coupling model. We have found that three coupling forms cannot lead 
to the well-behaved solution of the coincidence problem. Therefore, we introduce the forth 
coupling form in order to provide the ability to solve the coincidence problem. The properties 
for each coupling model are summarized as follows. 

For coupling (I), Q = y^2/3K(3p c (X + 3HX), the fixed points A are not stable. Their 
eigenvalues can be reduced to case of non-interaction. The effects of the interaction in the 
evolution of the energy density parameters are discussed in both qualitative and numerical 
analysis. We find the stability of the fixed points B numerically. By using three runaway 
potentials including exponential, Gaussian and inverse power law potentials, the results are 
that the fixed points B are stable and the stability is independent of sign of the coupling. 
However, they cannot solve the coincidence problem since they give the stable fixed points 
at which Ox = 1- However, it is instructive to study the effect of interaction in this coupling 
model. This effect is examined and the result is shown in Figure 1. For the exponential 
potential, V = Voe _??:r , the standard evolution of the universe, which has the matter- and 
radiation-dominated period, is obtained when y^r) > 0. It corresponds to the runaway poten- 
tial since the sign of y is the same as the sign of x. For the inverse power law and Gaussian 
potentials, both = 1 and y* = — 1 give the well-behaved evolution when ij > — 1 and 
T] > respectively. For the fixed points C, they can make an explanation of the coincidence 
problem. However, the coupling parameter (3 is of order 2, which makes the evolution at the 
early time different from the actual evolution as shown in the Figure 2(b). 

For coupling (II), Q = aHp c , the fixed points A are not stable. Their eigenvalues 
cannot be reduced to case of non-interacting one. This is because there is av 2 /2y 2 term 
which does not vanish when a = 0. For the fixed points B, the stability depend upon sign of 
the coupling. The singularity appears in the case of a > at the matter-dominated period, 
but is not in the case of a < for all potentials. This is because a > corresponds to energy 
transfer from dark matter to dark energy. This means that the energy density of dark matter 
increases when we go back in time. During matter-dominated period, the value of energy 
density of dark matter must take the enough value in order to give the observed value at the 
present time. For this reason, the singularity is emerged when v — > 1 and y — > 0. On the 
other hand, a < corresponds to energy transfer from dark energy to dark matter. There 
is no a singularity in this case because the value of the energy density of dark matter does 
not need to be large in the past. This implies that the energy transfer is allowed from the 
dark energy to dark matter only. Moreover, there are the other fixed point corresponding to 
a = — 3. However, a is too large in value for the evolution to agree with the cosmological 
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observations. Moreover, it gives non-dynamical evolution of the dark matter energy density. 

We also show that the suitable potential for the well-behaved evolution of three-form 
dark energy is power law potential. This property can be seen by comparing the evolution 
behavior of wtot for exponential potential (Figure 3) and inverse power law potential (Figure 
4). There is the oscillating behavior of wtot f° r exponential potential. We interpret this 
behavior from the last term in equation (3.20). The enormous value of x will factor this 
term to dominate even though z 2 approaches to zero. It is important to note that the inverse 
power law potential may give an unstable evolution of the cosmological perturbations since 
it provides negative sound speed square [20]. It is of interest to investigate whether the 
interaction can alleviate this behavior and we leave this investigation in further works. 

For coupling (III), Q = Tp c , One cannot eliminate the Hubble parameter H from the 
dynamical equations. We introduce new dimensionless variable in order to eliminate H by 
s = Hq/H. The results from this coupling are very similar to ones of the coupling (II), 
just a becoming 7s. The difference is that a is constant throughout the evolution of the 
universe, while 7s varies with time. However, by comparing to the coupling model (II), there 
is an advantage of this coupling model. This coupling model provides a tiny effect of the 
dynamical evolution in the past since s — > 0. This provides the wider range of the coupling 
parameter and leads to the very slight deviation of the evolution from the standard evolution 
of the universe in the past as shown in Figure 5. 

For coupling (IV), Q = V6T(X+3HX)/k, all of the fixed points A are not stable. Their 
eigenvalues can be reduced to the non-interacting case. However, we do not explicitly show 
the eigenvalues because of their complication. For the fixed points B, they can give the stable 
fixed points at which Qx/^c ~ 7/3 as shown in Figure 7. This provides the ability to solve 
the coincidence problem while the evolution is in agreement with the standard evolution of 
the universe. 

Even though there is no theoretical motivation for this coupling form, it may be in- 
terpreted as a physical interaction associated with particles since it does not depend on the 
expansion rate of the universe. T may also be interpreted as a decay rate like T in coupling 
model (III). This is one of the advantages for this coupling model. 

There is the cosmological constant fixed point, which A = 0, for coupling model (I), (II), 
(III) and non- interacting case but it does not exist for coupling model (IV). The important 
properties of this fixed point are that it is the stable fixed point with px = constant and 
ujx = constant = — 1. This is the property of cosmological constant dark energy and leads 
to the result that there are no stable fixed points at which Qx/^c ~ O(l). Covering the 
cosmological constant dark energy by setting A = or V = Vq is one of the important features 
for the three-form dark energy model which is distinguished from scalar field dark energy. 

By comparing the results to the interacting scalar field with type (II) and (III) [22, 23], it 
gives the similar result in which the coincidence problem will be solved. However, there is no 
negative energy density of the dark energy in our coupling model. We also note that there is 
only one parameter 7 to adjust in our coupling model while it has to have two parameters, cu 
and a c in the scalar interacting case. This suggests that it is easier to provide the possibility 
to solve the coincidence problem in interacting three-form field. One of the advantages of 
our model is that it automatically provides wtot slightly less than —1 before it reaches the 
stable fixed points as shown in Figure 7. This leads to the convenient way to fit the result 
with the observational data since recent observations suggest that uj = —1.1 ±0.14 [29]. 
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